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Abstract 

In time-resolved Brillouin scattering (also called picosecond ultrasonic interferometry), the 
time evolution of the spatial Fourier component of an optically excited acoustic strain distribution 
is monitored. The wave number is determined by the momentum conservation in photon-phonon 
interaction. For linear acoustic waves propagating in an homogeneous medium, the detected time- 
domain signal of the optical probe transient reflectivity shows a sinusoidal oscillation at a constant 
frequency known as the Brillouin frequency. This oscillation is a result of heterodyning the 
constant reflection from the sample surface with the Brillouin-scattered field. Here, we present an 
analytical theory for the nonlinear reshaping of a propagating, finite amplitude picosecond acoustic 
pulse which results in a time-dependence of the observed frequency. In particular, we examine the 
conditions under which this information can be used to study the time-evolution of the weak-shock 
front speed. Depending on the initial strain pulse parameters and the time interval of its nonlinear 
transformation, our theory predicts the detected frequency to either be monotonically decreasing or 
oscillating in time. We support these theoretical predictions by comparison with available 
experimental data. In general we find that picosecond ultrasonic interferometry of nonlinear 
acoustic pulses provides access to the nonlinear acoustic properties of a medium spanning most of 
the GHz frequency range. 
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1. Introduction 

In picosecond laser ultrasonics ultrashort optical pulses (typically of hundreds of 
femtoseconds duration) are used both to generate and to detect coherent acoustic pulses. 1 ' 2 In 
particular, the so-called, picosecond acoustic interferometry 1 ' 3 ' 4 is a powerful opto-acousto- 
optical method of nondestructive testing and evaluation of optically transparent materials in 
the GHz-range of acoustic frequencies. It consists in the generation of a picosecond acoustic 
pulse by a femtosecond laser pulse through laser action on light absorbing optoacoustic 
transducer and subsequent monitoring of the propagation of this acoustic pulse in the bulk of 
the transparent material by time-resolved reflectivity measurements. This is possible because 
of the interaction between the probe femtosecond laser pulse with picosecond acoustic pulses 
due to the acousto-optic effect. Governed dominantly by the law of the momentum 
conservation in the photons interaction with acoustic phonon, i.e., in the Brillouin scattering, 
picosecond ultrasonic interferometry also called time-resolved Brillouin scattering , monitors 
the time evolution of the spatial Fourier component of an optically excited acoustic strain 

distribution with wave number q = Ak . Here, Ak is the modification of the optical light 
wavevector in the probed medium. In the small-amplitude limit, a picosecond acoustic pulse 
propagates at a constant velocity, c 0 , without nonlinear changes to its shape. The detected 

time-domain signal of the optical probe transient reflectivity dR(t)IR, which includes a 
constant reflection from a sample surfaces/interfaces against which the Brillouin- scattered 
field is heterodyned (Fig. 1), shows a sinusoidal oscillation at a constant cyclic frequency 

a> dRIR equal to the Brillouin cyclic frequency <x> B = qc 0 . With Ak known, the measurement 

allows determination of the velocity c 0 (co B ) of acoustic phonons at frequency co B . 

Measuring sound velocity through the Brillouin frequency detection provides additional 
information necessary for the determination of the thickness of thin transparent layers in 
picosecond ultrasonics 8 ' 9 , where the detection either of the arrival times of the subsequent 
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acoustic echoes ' or the subsequent times of the abrupt amplitude and/or phase variation of the 
Brillouin signal 7 ' 10 ' 11 provides opportunity to determine the ratio of the film thickness to the 
sound velocity, but not separately the film thickness and the sound velocity. Time-resolved 
Brillouin scattering experiments have been extended for the measurements of the velocity of 
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shear hypersound . Performing picosecond ultrasonic interferometry at two different angles of 
probe incidence it is possible to determine simultaneously the sound velocity and the optical 
refractive index of the transparent material 13 ' 14 . Recently time-resolved Brillouin scattering 
experiments has been conducted in spatially inhomogeneous media 7 ' 15 . By measuring how the 
Brillouin frequency varies with the time delay of the probe laser pulse, i.e., with the propagation 
time of the acoustic pulse, it is possible to determine with nanometers scale spatial resolution the 
in-depth variation of sound velocity in the material ' .By performing measurements at few 
angles of probe incidence, it is possible to determine simultaneously the in-depth profiles of both 
sound velocity and of the optical refractive index 15 . Moreover from the measurements of the 
variations in time of the Brillouin signal amplitude it has been possible to determine the spatial 
profile of the distribution of the photoelastic modulus. 15 Recently developed theory 16 has 
predicted important advantages that would be achieved with this new technique of the nanoscale 
in-depth imaging if the amplitude and the phase of the transient optical reflectivity are monitored 
simultaneously using ultrafast time-resolved optical interferometry ' , which provides 
opportunity to separate precisely through measurements, and not by numerical signal processing, 
the time-varying amplitude and the time-varying frequency of the Brillouin signal. The 
perspectives of combining measurements with s-polarized and p-polarized laser probes has been 
also discussed. 16 Under the particular circumstances, for example for the depth-profiling of the 
shear rigidity of nanomaterials, it would be advantageous to monitor the Brillouin scattering 
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from the propagating nanoacoustic pulses by time-resolved ultrafast optical polarimetry. 
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Here we develop the theory of the time -resolved Brillouin scattering from the nonlinear 
coherent picosecond acoustic pulses, i.e., from the acoustic pulses of the finite and not of the 
infinitely small amplitude. In such pulses nonlinear elastic interactions between different 
components of the wide-frequency-band spectrum are important at the temporal/spatial scales of 
the experimental observation. These synchronous acoustic interactions lead to the nonlinear 
transformation of the picosecond acoustic pulse profile and the formation of the so-called weak 
shock front or several weak shock fronts, depending on the initial profile of the laser-generated 
picosecond acoustic pulse. We present an analytical theory describing how nonlinear reshaping 
of a propagating, finite amplitude picosecond acoustic pulse results in a time-dependence of the 
observed frequency a> dRIR {t) = co B (t) . In particular, we establish the conditions under which 

^dR/R^) = ( l v shi t ) where v sh {t) is time-evolving speed of the weak-shock front that can be 
formed through the nonlinear acoustic transformation of the picosecond strain pulse. Depending 
on the parameters of the initial picosecond strain pulse and the monitored time interval of its 
nonlinear transformation, our theory predicts the frequency co dRI R {t) = co B {t) to either be 

monotonically decreasing or oscillating in time. As the weak shock velocity depends on the 
nonlinear acoustic parameter of the material then the time-resolved Brillouin scattering provides 
access to the measurements of the acoustic nonlinearity in the GHz frequency range. The 
possible influence of the linear absorption of acoustic waves on the nonlinear Brillouin 
frequency is also analyzed. The theoretical predictions are compared with the available 
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experimental results for the bipolar and uni-polar nonlinear acoustic pulses. 

It is worth underlining, starting from here in the Introduction, that in general the 
technique of time -resolved Brillouin scattering monitors the dynamics of a particular spatial 
Fourier component of the acoustic field and not of a particular frequency of the acoustic field. 
Only in the case of the linear acoustic field the measurements provide information on the velocity 
of the acoustic phonon with the wave-number equal to this particular spatial Fourier component 
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(and with the frequency equal to this particular Brillouin frequency). In the case of the nonlinear 
acoustic fields it is not the case. The theory predicts that the detected frequency in the case of light 
scattering by nonlinear picosecond acoustic pulses is time-dependent and the experiments are 
providing information rather on the nonlinear phenomena in the complete frequency spectrum of 
the acoustic pulse than on the material elastic nonlinearity at a single time-dependent measured 
frequency. This is because the nonlinear variations in the pulse profile and, consequently, in its 
spatial Fourier spectrum are due to the interaction between all the spectral components of the 
acoustic pulse. The frequency of the detected oscillating component in transient reflectivity 
depends not only on the amplitude of this spectral component but also on the amplitudes and 
relative phases of the other components of the spectrum. The most obvious illustration of this 
statement is the case where the detected oscillation frequency is proportional to the velocity of the 
shock front, which is controlled by the amplitude of the acoustic pulse and not by the amplitude of 
a particular frequency component. In general, if the initial amplitude of the Brillouin frequency 
component is the same in two nonlinear acoustic pulses of the different profiles (of the different 
spectra) the variation in time of the measured frequency will be different. In other words, 
picosecond ultrasonic interferometry does not monitor the nonlinear self-action of a particular 
time-dependent frequency component, but nonlinear multi-frequency, multi-scattering acoustic 
processes in picosecond elastic pulses. Thus, picosecond ultrasonic interferometry of nonlinear 
acoustic pulses provides access to the nonlinear acoustic properties of a medium covering most of 
the GHz frequency range. 

The article is structured as follows. In Section 2 we remind the backgrounds of signal 
formation and description in picosecond ultrasonic interferometry (time-resolved Brillouin 
scattering). Section 3 is devoted to presentation of theory developed for the case when the acoustic 
signals, probed by picosecond ultrasonic interferometry, are nonlinear. In particular, in Section 3.1 
the transient optical reflectivity signals are predicted at the times scales preceding the formation of 



6 

weak shock front in a uni -polar nonlinear acoustic pulse. In Section 3.2 the manifestations in the 
transient reflectivity signals of the formation and propagation of infinitely thin shock fronts are 
described (regime of high acoustic Reynolds numbers). The theory presented in Section 3.3 
accounts for the finite duration of the weak shock fronts caused by the frequency-dependent 
absorption of the hypersound. In Section 3.4 a particular case of transient broadening of weak 
shock fronts, initially generated by ultrafast laser action, is analyzed. Section 3.5 is devoted to the 
analytical description of the Brillouin signals in a most general case of multi-polar acoustic pulses 
and arbitrary acoustic Reynolds numbers. Finally in Section 3.6 an important particular case of 
symmetric bipolar acoustic pulses is studied in more details. Section 4 presents two applications of 
the developed theory to the analysis of the available experimental results. In Section 4.1 we 
interpret earlier reported experimental results obtained by spatially resolved Brillouin scattering 
technique. In Section 4.2 we use the derived theoretical formulas for the analysis of the available 
experimental results obtained by time-resolved Brillouin scattering technique. The conclusions are 
summarized in Section 5. 

2. Picosecond ultrasonic interferometry of coherent acoustic pulses 

In picosecond laser ultrasonics ultrashort probe laser pulses are scattered by laser- generated 
picosecond acoustic pulses due to the acousto-optic effect. 6 ' 26 ' 27 Thus, by monitoring the reflected 
light, the information on the picosecond acoustic pulses can be obtained. 1 In the most of the 
experiments the transient relative changes dR(t) I R in the intensity of the reflected probe light are 

detected, where R is the sample reflectivity in the absence of pump laser action. The optical 
intensity reflectivity R and electromagnetic field reflectivity r are related by 



dR(t) / R = 2Re[dr(t) / r] . 



(1) 
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The contribution to this signal, which is caused by acousto-optic interaction of plane 
collinearly propagating acoustic and light waves in isotropic medium, is proportional to the photo- 
elastic (acousto-optic) coefficient and the overlap integral of the acoustic strain field with the 
optical field sensitivity function. 1 The integral is over complete volume probed by the light and the 
sensitivity function depends on the specific geometry of light scattering in the sample in the 
absence of the laser-generated coherent acoustic waves. 1 ' 5 ' 10 In particular, when the light normally 
incident on the sample probes the propagation of the acoustic pulse in optically homogeneous half 
space (Fig. 1) the photo-elastic contribution to Eq. (1) is given by 



r 



= -ic[j 1zz (t,z)e i2b dz, (2) 



where C is the known complex constant 1 ' 5 ' 10 , Tj zz (t,z) is the strain field in the longitudinal 
acoustic wave propagating along the z axis and k is the complex wave number of the probe light in 
the medium. In accordance with Eq. (2) the time-resolved measurement of the reflectivity provides 
information on the evolution in time of the Laplace component of the acoustic field with the 
complex-valued Laplace parameter p = -ilk . 5 The solution in Eq. (2) is particularly simple, when 
the strain field is composed of the coherent waves propagating in one direction. For example, the 
strain pulse propagating in the positive direction of z axis can be presented in the form 



Tj zz (t,z) = T] zz (t,^ = z-c 0 t) = —\ ff z ,(t,q)e^dq. (3) 

2k 



Here c 0 is the velocity of the linear acoustic waves, i.e., the velocity of the acoustic waves of 
the infinitely small amplitude, and fj zz {t,q) describes slow evolution with time of the spatial 
Fourier spectrum of the coherent acoustic pulse. Note that J; = z-c Q t is the moving spatial co- 
ordinate, and, so, the slow time t in rj zz (t,^) and in rj zz (t,q) describes the evolution of the strain 
field in the co-ordinate frame moving with linear sound velocity. When the strain pulse (3) is 
completely inside the probe half-space z > 0 , i.e., when the value r] z Xt, z = 0) predicted by Eq. (3) 
is negligible, the integration in Eq. (2) can be extended to complete space 



dr 
r 



= -iC\ + j zz (t,z)e i2kz dz. (4) 
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In the optically transparent medium the value of k is real, and, thus, Eq. (4) predicts that the 
time-resolved measurements of the reflectivity monitor the evolution in time of the component of 
spatial Fourier spectrum of the acoustic wave with the wave number q equal twice the optical wave 
number. The substitution of Eq. (3) into Eq. (4) results in 



rlr 

— = -iCri zz (t,q = 2k)e i2kc «. (5) 
r 



The relation q = 2k is a well-known result of the momentum conservation in back-scattering 
of photon by phonon ' . Thus, mathematically, it is the spatial Fourier transform that ensures the 
momentum conservation in the photons-phonon interaction. Because of the relation q = col c 0 , 

where co is the cyclic frequency, it is common to say that time resolved reflectivity measurements 
monitor the time dynamics of the spectral Fourier component of the coherent acoustic field at 
frequency co = 2kc 0 = co B , which is known as the Brillouin frequency or the Brillouin frequency 
shift of the light scattered by acoustic phonon. The measurement of the Brillouin frequency 
provides opportunity to find the acoustic velocity in the material through the relation c 0 = co B /(2k) . 
These statements can be straightforward verified by presenting the strain pulse in the form 

I r +co 

V zz {^z)^Ti u (t,T = t-zlc 0 ) = — \ ij u (t,oS)e WJT dco, (6) 

2n J -° 0 



where r = t-z/c Q is the so-called retarded time 5 and r) u (t,co) is the slowly varying in time 
Fourier frequency spectrum of the acoustic pulse. Substitution of Eq. (6) into Eq. (4) results in 



dr 

— = -iCfj zz (t,co = co B )e ia "' t . (7) 
r 



The solution in Eq. (7), which is equivalent to one in Eq. (5), confirms that time-resolved 
measurements of the probe light intensity reflectivity (see Eq. (1)) provide access to the detection 
of coherent acoustic waves at the Brillouin frequency. This is the reason to call this domain of 
picosecond laser ultrasonics not only "picosecond ultrasonic interferometry" 1 but also sometimes 
to call it "time-resolved Brillouin scattering" . However the fact, that the Brillouin frequency can 
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be measured using time -resolved variations of the optical reflectivity caused by laser-generated 
coherent acoustic fields, does not mean that, by measuring the frequency of the oscillations in the 
dR(t)/R signal, we get exactly the Brillouin frequency <x> B = 2kc 0 . By separating the amplitude 

and the phase in the complex C and rj zz (t, co) , i.e., C = A c exp(i(p c ) , rj zz (t, a> B ) = A^(t)exp^tp (t)\, 

and substituting Eq. (7) into Eq. (1), we derive 

dR(t) I R = 2A c A t] sin(co B t + <p ri +<p c ). (8) 

From Eq. (8) it follows that the transient reflectivity signal induced by photoelasticity 
oscillates at frequency 

*W = 5 (<*y + <p n +<p c )/dt = co B + dq> n I dt , (9) 

which is, in the general case, different from the Brillouin frequency defined above, i.e., 

Below we will evaluate co dR/R in the particular case where the time -resolved Brillouin 
scattering is applied for the monitoring the propagation of the nonlinear acoustic pulses. We will 
describe what information can be extracted from the measurements of co dRIR . In particular, it will 
be demonstrated that under some conditions 

co dR/R /(2k) = (co B + d<p v I dt)/(2k) = c 0 + (l/(2^))5^ / dt provides the information on the velocity 

of the weak shock front in the profile of the coherent acoustic pulse. It will be also discussed how 
the parameter of the acoustic nonlinearity of the medium could be accessed through time-resolved 
Brillouin experiments. 

Below we will show that dqt^ I dt could be non-zero and could provide non-zero 

contribution to the experimentally measured co dR/R because of the elastic nonlinearity of the tested 

medium. So it is tentative to call co dR/R = co B + d<p v I dt just the time-dependent nonlinear Brillouin 

frequency co dRIR = co B + dq> n I dt = co B (t) (or nonlinear frequency shift). In principle this 

terminology is acceptable. However, when using this terminology it should not be forgotten that 
linear Brillouin frequency co B is the characteristic of the tested material, which does not depend on 
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the parameters of the monitored picosecond acoustic pulse, while, as it will be demonstrated below, 
d(p v I dt and consequently also co dRIR depend additionally on the temporal profile, i.e. on the 

complete frequency spectrum, of the monitored picosecond strain pulse. 



3. Picosecond ultrasonic interferometry of nonlinear acoustic pulses 



Propagation of plane longitudinal acoustic pulses of finite amplitude is controlled in the 



frame of quadratic approximation of the nonlinear acoustics 23 by the Burgers equation 



do do b d 2 o . 

— + so r = 0. (10) 

dt d% 2p Q dC 



Here o(t, E, = z - c Q t) describes slow evolution in time of the profile of the particle velocity, 

which is evaluated in the sonic coordinate frame. This evolution takes place due to the effects of 
quadratic nonlinearity, which are accounted by the second term in Eq. (10), where a is the 
nonlinear acoustic parameter , and due to sound absorption proportional to square of frequency, 
which is accounted by the third term in Eq. (10), where b is the effective viscosity 5 and p 0 is the 
equilibrium density. The dispersion of the acoustic velocity is neglected in Eq. (10). 

Let us first present a simple particular example demonstrating the inequality co dRIR ^ co B 

for the nonlinear acoustic pulses described by Eq. (10). We note preliminary that, because particle 
velocity and the longitudinal strain are related to particle displacement u in the acoustic wave by 
o = dul dt and Tj a =du/ dz , respectively, their relation in the moving coordinate frame is 

O = -C 0 7J a .So 

fj zz (t,q = 2k) = --f v(t,Z)e i2k *dt. (11) 

/ • J— CO 

c 0 



This relation will be used in the following for the evaluation of the amplitude and the phase 
of the complex strain spectrum fj zz (t,q = 2k) entering the solution for transient reflectivity in Eq. 
(8) and Eq. (9). Without loss of generality it can be considered that the acoustic pulse is located 
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around % = 0, i.e., around z = c 0 t . Then for the short acoustic pulses with the spatial length l a (t) , 
satisfying the inequality 2kl a (t)«\, i.e., l a {t) «A.I{Att) , where X is the optical wavelength in 
the medium, Eq. (11) can be approximated by 



fj a {t,q = 2k)« v(t,&(l + 2ikt)dZ * eap/— f v(t,&& 



(12) 



poo 

where S = v(t, <^]d<^ denotes the "surface area" of the acoustic pulse velocity profile, 

J— oo 



23 



which is known to be constant in nonlinear acoustics due to acoustic momentum conservation. 
The obtained relation (12) demonstrates that, as far as the profile of the acoustic pulse changes in 
time, i.e., du/dt^O, the phase of the complex strain could contain a time-dependent part 



2k 



poo 

u(t,^)^d^ and co dRIR * co B could hold. Moreover, using the Burgers equation (10) and 

J— oo 



S J 

applying the integration by parts, the deviation of co dRIR from co B can be reduced to the form 



03 dRIR ~ m B 



^ 2k r° dv(t,£) 2k p 

~ S J- dt S > 



do b d v 
oo h 



d% 2 Po dC 7 



^ = ^[u 2 d{. (13) 

C J-co 



In the derivation of the final result in Eq. (13) we have used the condition of the absence of 
the acoustic field at £ = ±qo . The solution in Eq. (13) demonstrates that the deviation of the 

measured frequency co dRIR from the Brillouin frequency co B is caused only by the acoustic 
nonlinearity, i.e., is absent in the absence of nonlinear processes (in case s = 0). In the absence of 
nonlinearity linear absorption does not induce the difference between co dRIR and co B . In the 

particular case of short acoustic pulses co dRIR - a> B is proportional to the ratio acoustic pulse energy 

pco poo 

and the acoustic pulse momentum co dRIR - co B ~ sk\ v d% I vd% . This relation indicates that 

' J— CO J— GO 

co dRIR -co B could contain information on the product su a of the parameter s of the acoustic 
nonlinearity and the characteristic amplitude v a of the particle velocity in the strain pulse. 

In the general case of the acoustic pulse characterized by the spatial length l a , i.e., the 
duration z a =l a /c 0 , and the amplitude v a , the evolution of the pulse profile predicted by Eq. (19) 
is known to be controlled by the parameter s Re , where the notation Re is introduced for the so- 
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called acoustic Reynolds number 23 Re = 2p Q c Q o a T a lb . If the introduced parameter is small, 

s Re « 1 , then the dissipation dominates over the nonlinearity, and the nonlinear corrections to 
the acoustic pulse profile are small and can be taken into account by the method of the successive 
approximations starting from the solution of the linearized Eq. (10). Moreover if the parameter 
e Re is small in the initial instance of time t = 0 , then it will be small later, and the nonlinear 
effects will be small over complete observation time . This simple regime of acoustic pulse 
propagation is not of prime interest of our investigation and is not discussed below in details. We 
are interested here in the analysis of the situation, when initially the nonlinear effects dominate, i.e., 
e Re(£ = 0) » 1 , and in the first approximation the evolution of the pulse profile can be described 
neglecting in Eq. (10) the absorption, i.e. by the so-called simple-wave equation 23 



— + so — = 0. (14) 

dt d% 



It is known, however, that in this case the nonlinearity leads to the formation of the so-called 
weak shock fronts in the profile of the acoustic pulse 23 . The duration of the shock front r sh (t) can 

be much shorter than the initial duration of the acoustic pulse. It means that the profile of the 
acoustic pulse can not be characterized by a single characteristic time. In the process of the shock 
steepening the local acoustic Reynolds number in the vicinity of the shock, Re = 2p 0 c 0 u a r sh (t) I b 

progressively diminishes. So this is the vicinity of the shock front, where the corrections to the 
pulse profile due to the dissipation should be first taken into account in the regime characterized by 
s Re(t = 0) » 1 . Thus the common way to describe the pulse evolution in the regime 
sRe(t = 0) »1 of interest here is first to analyze the pulse evolution, predicted by Eq. (14) and 

23 

then to take into account the influence of sound absorption , which is the most pronounced in the 
vicinity of the weak shock fronts. 



3.1. Probing transformation of nonlinear acoustic pulse in ideal media before weak 
shock front formation 

In ideal media, characterized by negligible effective viscosity the propagation of the 
nonlinear acoustic pulses is described by the simple-wave equation (14). The transformation of the 
pulse profile, predicted by the simple- wave equation (14) is easy to "visualize". 28 " 30 In the 
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coordinate system (l>,£) all the points (t> (0) = u(0, ^, (0)), ^ (0)) of the initial profile v(t = 0, £,) 
are moving with increasing time horizontally following the law 
(Uj(0 = u,.(0),£(0 = £(0) + 6U,(0)0 . 28 ~ 30 In Fig. 2 (a) we illustrate how this motion of the profile 
points leads to the formation of the shock front in the initially triangular pulse. The lengths of the 
horizontal arrows in Fig. 2 (a) increase as £L>,(0)£ . The letters A, A' and A" denote in Fig. 2 (a) 
the positions of the acoustic pulse peak in the initial time moment t = 0 , in the intermediate time 
moment 0 < t < t sh before the shock front formation and in the moment t = t sh of the front 
formation, respectively. If the initial full-width-half-maximum (FWHM) length of the triangular 
pulse is l a and the amplitude is v a , then the time for the shock formation is t sh = l a l(su a ) . This 
formula can be used for the estimation of the shock formation time in the pulses of more 
complicated shape if their profiles can be characterized by only two scales, i.e., l a and v a . 
Moreover the transformation of the pulse profile presented in Fig. 2 (a) captures the most essential 
features of the transformation for 0 < t < t sh of the single-polar pulses of other shapes 
characterized by only two scales. The analytical description of the pulse profiles in Fig. 2 (a) is 



0, if $<-2l a \ 

u a (2l a +4)/(l a +eu a t), if -2l a <4<-l a +euj; 
uJ/(-l a +su a t), if -l a +£UJ<Z<0; 
0, if 0<4. 



The substitution of Eq. (15) into Eq. (11) results in 



^ (t a = 2k) = - f sin N a -£*V)] e *. _ sin[fc(* fl +ay)] ^ ) _ 2 ^ +;fe(V 
2ikc Q [ k(l a -so J) k(l a + so J) J 



For the picosecond laser ultrasonics the most typical and the most important is the situation 
where the initial length of the photo-generated acoustic pulse is much shorter than the optical 
wavelength. In this case the first approximation to the precise solution in Eq. (16), which includes 
the terms of the first order in small parameter 2kl a « 1 , is 
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r/ zz (t,q = 2k) = -^ L (l-2ikl a +2iksvj/3) = -^e " 3 " =-^- e 3 Wj .(17) 



2 , ■, , 2iH, 

V I V I -2ikl a +i-ksv a t {J I - 2 <*'«+- 



With the accepted precision, the derived solution in Eq. (17) predicts that at the stage of the 
shock front formation in the profile of the picosecond acoustic pulse the amplitude of the 
oscillation in transient reflectivity is a constant independent of the material nonlinearity, while the 
deviation of the oscillation frequency from the Brillouiin frequency is a constant proportional to 
acoustic nonlinearity e , acoustic Mach number M a =v a l c a and the Brillouin frequency itself 



2k 1 

G>dRIR-°>B~—tVa=^ £M a 0 >B- (18) 



It can be verified that, for the assumed short triangular acoustic pulse, the result in Eq. (18) 
can be much easier obtained using Eq. (13). 



3.2. Picosecond ultrasonic interferometry of acoustic pulse containing weak shock front 

In Fig. 2 (b), by the same method of geometrical transformation, we illustrate that after 
the shock formation, i.e., at times t > t sh the profile of the acoustic pulse predicted by Eq. (14) 
becomes multi-valued, i.e., non-functional. To avoid this unphysical prediction, the shock front 
should be introduced in such position % = £ sh (t) that the surface area of the pulse is conserved . 

This geometrical construction predicts that for t > t sh the length of the strain pulse starts to 
increase, l a {t) = l a + g sh (t)/2 , while the amplitude of the acoustic pulse starts to diminish 
v a it) = uj a I l a it) . The letters A" and A"' denote in Fig. 2 (b) the positions of the acoustic pulse 
peak in the moment t = t sh of the front formation and in some moment t > t sh after the shock front 
formation, respectively. The analytical description of the pulse profile takes the form 
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where 



u.(0 = «./Jl + 



sM, 



t-t. 



= vl 



''a J 



1 + 



r t-t,^ 

sn 

2t 



(20) 



1 + 



'a J 



= 21 



1 + 



2f 



(21) 



//(£) is the unit step function ( H(i; < 0) = 0, H(i; > 0) = 1 ) and sgn(£) is the sign function 
( sgn(<^ < 0) = -1, sgn(^ > 0) — 1). It can be checked that the velocity of the shock front propagation 
in the moving coordinate frame is controlled by the nonlinear acoustic parameter and the 
amplitude of the acoustic pulse, v sh (t) = d% sh (t)l dt = su a (t)l2 . 
The substitution of Eq. (19) into Eq. (11) results in 



fj (t,q = 2k) = 



2ikCr, 



,2*&»(f) 



\ \l _ e -4*['„+^(')/2]ll 

4ife[/ fl +^(0/2] L f 



(22) 



At the time scale, when the acoustic pulse, even broadened by the nonlinear processes, is still 
importantly shorter than the optical wavelength, i.e., 2kl a (t) = 2k[l a + % sh {t) 1 2] « 1, the solution 
in Eq. (22) can be approximated by the formula 



rj zz (t,g = 2k) = -^e 3 3 , 
c 0 



(23) 
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which predicts that after the shock formation the difference between the detected frequency and 
the Brillouin frequency starts to diminish 



CO 



4 

3 



dRiR -®b~- kv, h (t) = - sM a (t)co B 



1 

3 



(24) 



In 



Eq. (24) M a (t) denotes 



M a (t) = M a /ll + 



t-t 



sh 



V T a J 



= M I Jl + 



the current acoustic Mach number, 
The condition 2kl a it) « 1 , required for 



the validity of Eq. (23) and Eq. (24), can be rewritten under the assumption 2kl a (t =0) «1 as 
t — t sh « 2\2kl a \ l t sh = . Because of t { ^»t sh it can be concluded that the considered 
approximation could be valid at time scale significantly exceeding the characteristic nonlinear time 
t sh of the shock front formation. 

With the increasing time, due to the nonlinear broadening of the acoustic pulse, the role 
of the second term in the figure brackets in Eq. (22), which describes the light scattering by the 
trailing edge of the strain pulse (see Fig. 2 (b)), progressively diminishes. This term is comparable 
with the first term only if 2kl a {t)<\. At larger times the solution in Eq. (22) predicts that the 

dominant scattering of the probe light is by the leading shock front of the acoustic pulse, where the 
strain gradients are much higher than in the trailing front 



rj (t,q = 2k) = - ^ e 2ik ^ (t) = ^ e 2 *** (0+iW2 



2ikc n 



2kc n 



(25) 



In this regime the amplitude of the oscillation in transient reflectivity diminishes with time 
following the fall in the pulse amplitude, described by Eq. (20), while the oscillation frequency is a 
measure of the shock front velocity. In fact, it follows from Eq. (25) that 
co dRIR -co B « 2kv sh (t) = (sM a (t) 1 2)co B . So, if to the experimental measurement of frequency we 

associate the measurement of the velocity through the relation a> dRIR = 2kc dRIR similar to 

<x> B = 2kc 0 , which holds for the Brillouin frequency, then in the considered regime the measured 

velocity c dRIR is c dRIR = c Q + v sh (t) = c sh . Thus we are measuring the velocity c sh of the weak 



17 



shock front in laboratory, i.e., immobile, coordinate frame. The velocity c sh of the compression 

shock, i.e., of the one with u a (t)>0,rj zz < 0 , in medium with positive parameter of acoustic 

nonlinearity ( s > 0 ), is supersonic ( c sh > c 0 ). 

The regime of monitoring of the shock front velocity by time-resolved Brillouin 
scattering can be understood more generally, by transforming Eq. (11) through the integration by 
parts, into 



The solution in Eq. (26) has a clear physical sense that scattering of light is caused by strain 
gradients and not by the homogeneous strain. Homogeneous strain distribution can not contain the 
spatial Fourier component at q = 2k ^ 0 . When the nonlinear acoustic processes progressively 
shorten one front of the strain pulse and broaden the other one, high spatial Fourier frequencies are 
progressively more and more spatially localized in the narrower front. Finally, when the broader 
front of the acoustic pulse provides negligible contribution to Eq. (26), because of low spatial 
gradients of the strain in this front, then time-resolved reflectivity measurements monitor the 
evolution of the shorter front, i.e., of the shock front, only. The formal substitution in Eq. (26) of 
the delta-localized shock front, i.e., du(t, « -v a {f)8{% — % sh {t)) , immediately reproduces 

differently derived solution in Eq. (25). 

We can conclude, that time -resolved Brillouin scattering always selects from the acoustic 
wave profile the spatial Fourier component with q-2k, however, in the case of the pronounced 
nonlinear transformation of the coherent acoustic pulse profile accompanied by the shock front 
formation, this particular component (Brillouin component) is spatially localized in the shock front 
and moves due to nonlinear processes not with the velocity of the linear acoustic waves c 0 but 

with the velocity of the shock front c sh = c 0 + v sh (t) . Correspondingly, the measured frequency of 

the oscillation in transient reflectivity differs from the Brillouin frequency that should be observed 
in the linear regime of the coherent acoustic pulse propagation. 

It is worth repeating here that the formation of the weak shock front and its subsequent decay 
in amplitude and deceleration in propagation velocity are all the results of nonlinear multi- 
scattering acoustic processes involving all the spectral components of the wide-frequency-band 
picosecond acoustic pulse. They all depend on the particular initial profile of the laser-generated 




(26) 
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picosecond strain pulse. The deviation of oo dRIR from co B is not just the result of the nonlinear self- 
action of the spectral component with the frequency equal to co B . It is well known in nonlinear 

23 

acoustics that in the initially sinusoidal acoustic wave of finite amplitude, in the case of elastic 
quadratic acoustic nonlinearity accounted by Eqs. (10) and (14), the weak shock fronts are formed 
inside each spatial period of the wave, but they all propagate at the velocity of the linear sound c 0 , 
i.e., they are immobile in the moving sonic coordinate frame. So if a purely sinusoidal strain wave 
at frequency a> B could be laser-generated (for example, in the process of the stimulated Brillouin 
scattering of the pump laser radiation 28 ) its nonlinear transformation in the media with quadratic 
elastic nonlinearity would not result in the shift of oo dRIR relative to co B . 

The solution in Eq. (22) describes the transition from monitoring of the both fronts of the 
acoustic pulse (Eq. (23)) to monitoring of the shock front only (Eq. (24)). Introducing the compact 
notations 4k[l a +g sh (t)/2] = 4kl a (t) = %7d a {f)l X = A(t) for the normalized length of the strain 
pulse and sin c(x) = sin(x) / x for the sine function, we rewrite the solution in Eq. (22) as 



(t, q = 2k) = - Jl - 2 sin c(A) + sin c 2 (A / 2)e 



2i'£f rfl (0-!arctan 



l-sinc(A) 
(A/2)sinc 2 (A/2) 



2kc, 



(27) 



The square root in Eq. (27), E(f) = -Jl - 2 sin c(A) + sin c 2 (A / 2) , describes transitional 
modulation of the Brillouin signal amplitude (Fig. 3 (a)), which is superimposed on the asymptotic 
dependence (t) <x v a (t) . The dependence on A(0 of the phase of the signal in Eq. (27) results 

in transitional modulation of the frequency of the transient reflectivity oscillation, proportionally to 

1 - sin c(A) 
l-2sinc(A) + sinc z (A/2) 

dependence a> dRIR <xv sh (t)a> B oc sM a (t)a> B 



0(A) = - — — — — : — 2 (Fig. 3 (b)), which is superimposed on the asymptotic 



CO. 



dRIR 



1- sinc(A) 
l-2sinc(A) + sinc 2 (A/2) 



sM a (t) 



1 - sin c( A) 
l-2sinc(A) + sinc 2 (A/2) 



co B 
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(28) 

The transitional modulations are due to the interference between the light waves scattered by 
the leading and trailing edges of the strain pulse. They disappear when the scattering of light by 
one of the fronts becomes negligible. In accordance with the analysis presented above the 
transitional regime described by Eq. (27) takes place when the inequality 2kl a (t) < 1 holds. This 

condition can be rewritten as t — t sh < 2[2kl a Y 2 t sh = , demonstrating that, because of the strong 

inequality 2kl a « 1 , the monitoring of the motion of the shock front only could be possible only 

after time delay which is importantly longer than the characteristic nonlinear time t sh . 



3.3. Influence on the time-resolved Brillouin scattering signals of finite shock front 
duration. 

Infinitely short shock fronts of the nonlinear acoustic pulses in Fig. 2 (a) and (b) are just 
the first approximation, which is provided by the precise solution of the simple-wave equation (Eq. 
(14)) to the solutions of the Burgers equation (Eq. (10)) in the case, where the initial pulse is 
characterized by the inequality s Re(7 = 0) » 1 . In this regime in the next approximation the finite 
duration of the shock front could be taken into account by replacing in the solution (19) the 



multiplier sgn(£ - % sh (t)) by tanh 



, where tanh denotes the hyperbolic tangent 



function and the width of the shock front A£ sh (t) depends on the shock amplitude 



A£,(0 = c 0 r sh = 2b/[tp 0 u a (t)]= [4/(*Re)]f a 



1 + 



V 2t sh j 



P It can be verified that, with thus 



defined shock front width, the local acoustic Reynolds number in the shock front region 
Re = 2p 0 c Q u a r sh (t)/b is equal to Ms . Consequently inside the shock the estimate <?Reocl is 

valid and the processes of the linear high frequency absorption causing front broadening are as 
important as the nonlinear processes leading to front shortening. Front broadens following the 
diminishing of its amplitude, A^ sh (t) oc \lv a {t) . The solution, which is obtained in this 

approximation, seems to be accurate enough as far as the shock front retains its physical sense, i.e. 
as far as it is sufficiently shorter than the total length of the acoustic pulse, 
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A^(0 = 



^Re 



L.1 + 



2t 



1 + 



2t 



(29) 



In the regime <?Re»l this inequality reasonably holds until complete absorption of the 
acoustic pulse. However it is known 23 that thus derived solutions can not be used till t = <x> , 
because in the adopted approximation the next order effect, consisting in the diminishing of the 
pulse amplitude v a (t) , when the shock front of the finite duration replaces the infinitely short 
shock front, is not taken into account. This, additional to one described by Eq. (20), diminishing of 
v a (t), leads to additional broadening of the shock front due to the relation Ag sh (t) oc \l v a {t) , 

forming the inverse loop. The described slowly accumulating effect is known to lead at long time 
scale to the symmetrical single polar pulses with equal duration of the leading and the trailing 
fronts, as we will demonstrate a bit later through the analysis of a particular precise solution of the 
Burgers equation. Here we note that the conditions of the confident validity of the method of the 
shock front installation can be found, by requesting that the installation of the shock front 
A£ sh (t) oc 1/ v a {t) in the pulse profile leads to the negligible changes in the pulse amplitude even if 

the increase in the pulse width l a {t) is neglected, i.e., if the front of the increasing width is 

installed in the fixed initial pulse profile. This transforms the validity condition in Eq. (29) into 
more restrictive one 



A^(0 = 



U1 + 



r t-t A 

1 l sh 



2t 



«l.<fl) = l., 



(30) 



which even for ^Re»l provides the upper time limit for validity of the solutions 
t - t sh « 2{s Re/ 4) 2 t sh = . We conclude that, because of the inequality e Re » 1 , the obtained 
solutions are valid at the time scale importantly larger than the characteristic nonlinear scale t sh of 
the pulse broadening due to nonlinearity. If the inequality (2kl a Y 2 < (^Re/4) 2 holds, then 
*sh < *sh an d m the time interval < t < t^f the influence of the front broadening on the transient 
reflectivity signal takes place when the influence of the light scattered by the trailing front of the 
pulse is negligible. 
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To capture the essential influence of the shock front broadening on the transient 
reflectivity signal, but avoiding cumbersome formulas, we approximate the tanh-type shock front 
by the linear shock front 



tanh 



1, if £<£ ft (0-A£ ft (0; 

-<r/A£ A (0 if ^(0-A^(0<^<^(0 + A^(0; 

-1, if Z>Z sh (t) + AZ sh (t). 



Correspondingly, the spatial derivative of the pulse profile inside the shock front, i.e., when 
^(0-A^(0<^<^(0 + A^(0, is given by du(t, &/d£ - ~[u a (t) /(2A^ sh (t)] . Then the 
evaluation of Eq. (26) leads to 



rj zz (t,q = 2k) = 



sin(2*A£ A ) 



2*^(0 



(31) 



The solution in Eq. (31) indicates that continuous broadening of the shock front is the source 
of the modulation of the transient reflectivity amplitude, which is superimposed on the asymptotic 
dependence A^it) cc v a {t) . Note, that because of the inequality (29), the solution in Eq. (31) 

indicates that this modulation is activated only after the disappearance of the modulation caused by 
the light reflection from the long trailing edge of the nonlinear acoustic pulse. It is important that, 
in the frame of the accepted approximations, the absorption causes shock front broadening but 
does not induce shock front displacement relative to % = £ sh (t), which was predicted earlier in the 
absence of absorption. This results in the absence of the modification of the phase of the signal in 
Eq. (31) except at the time instances defined by the solutions of 2kAg sh (t) = m , where n = 1, 2, ... 

is the integer number. At these time moments the phase of the Brillouin oscillation abruptly 
changes by n . If the phase is appropriately unrapped, then these changes do not result in any 
modulation of the oscillation frequency additional to the asymptotic dependence 
co dR/R -co B cc v sh (t)co B oc sM a {t)co B revealed earlier. 



3.4. Influence on the signals of transient broadening of the initially laser-generated 
weak shock fronts 
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It is important to note that the conditions of the validity of above derived description of 
the time-resolved light scattering by the coherent acoustic pulses could be broader than given by 
the inequality e Re(t = 0) » 1 . In fact, the description of the finite width of the shock front 



A£*(0 - c 0 r sh = 2b/[tp 0 u a (t)] = [4/(*Re)K 



1 + 



v ^sh J 



, accepted above, is suggested by the 



stationary shock fronts predicted by existing analytical solutions of the Burgers equation. 
Stationary shock fronts are the time-independent smoothed step-like wave profiles moving with 
constant velocity, different from the velocity c 0 of the linear sound. 23 So the above proposed 

solutions had not taken into account that in the transient problems of acoustic pulse propagation it 
could take finite time to smooth the initially step-like profile. Let us illustrate this statement by a 
typical example from picosecond laser ultrasonics. In picosecond laser ultrasonics in many cases 
the initially photo-generated, i.e., just after the thermalization in the material of the energy 
absorbed from the femtosecond laser pulse, single polar acoustic strain pulse has a highly 
asymmetric shape which is much closer to that described by Eq. (15) and in Fig. 2 at t = t sh than at 
t = 0 . In other words the initial pulse and its evolution should be rather described not by Eq. (15) 
and Eq. (19) but just by Eq. (19) with the shift of the initial time moment, i.e., t — t sh =>t. In this 

description the characteristic width of the initial pulse is approximately equal to the initial length 
of trailing front. In the case of the emission of the acoustic pulses in the transparent medium from 
the light absorbing medium this characteristic length l a could be controlled by the length of optical 

energy transfer to the material atoms in the opaque medium or by the thickness of the opaque 
medium, for example. The length of the leading front could be controlled by the distance of 
acoustic wave propagation either during the femtosecond laser action or during the time of the 
absorbed energy transfer to the material atoms, or by the characteristic roughness of the interface 
between opaque and transparent media, for example. The initial width of the leading front is 
negligible in comparison with the width of the trailing front. So we can try to treat it as shock front 
of the initial amplitude v a . However, it would be incorrect to replace the infinitely thin shock front 
in Eq. (19) by the front of the finite width 



A£ sh (t) = c Q r sh = 2b/[qo 0 u a (t)] = [4/(e Re)]/ a ^/l + 1 l(2t sh ) immediately, that is starting from t = 0 . 
This follows from the precise solution of the linearized Burgers equation, i.e., of the equation 

dv b d 2 u 



8t 2p Q dC 



= 0, with the initial condition v(t = 0, = v a [\ -sgn(£)] . This precise solution 
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reads u(t = 0, £) = v a [l - erf (£ I pbt I p 0 ) J and predicts that the broadening of the front width 

A£ abs caused by the high frequency absorption in the absence of the nonlinear processes follows 

the diffusion law A£ ahs (t) = ^2bt I p 0 . This could have been expected because the linearized 

Burgers equation is the linear diffusion-type equation. We can conclude that, if the nonlinearity is 
taken into account, the duration of the shock front could be estimated by the earlier proposed 
A% sh (t) only at the time scales when A% sh (t) < A% abs {t) . This inequality reduces to 

t>2t sh /( £ Rd2-l) = t x . (32) 

When <?Re » 1, the inequality (32) is satisfied starting from the critical time t x « t sh . This 
confirms the validity of the solutions, that were derived above (for <?Re »1 at the time scales 
t sh <t < 2(£ , Re/ 4) 2 £ s/i of interest, when the shock front displacement in the sonic coordinate frame 
could be appreciated) through the installation of the stationary front of the width A% sh {t) < A% ahs (t) . 
Note, that for the acoustic pulses containing initially the shock front, the characteristic time 
t sh =l a l(su a ) should be understood not as the time of shock front formation, but just as the 
characteristic time of the possible broadening of the initial acoustic pulse due to the nonlinearity. 
In accordance with l a (t) = l a ^]l + t/(2t sh ) at t = 6t sh the length of the initial triangular pulse could 
be doubled. 

When the inequality in Eq. (32) is not satisfied, then the duration of the shock front to be 
installed in the solution Eq. (19) should be A£; ahs (t) = ^2bt I p 0 , and the condition of the negligible 
influence of the shock front broadening on the front amplitude results in the equality 
■yj2bt I p Q « l a . This inequality can hold even for the acoustic waves of low amplitude, when 

t « t sh in the experimental observation time window, if t « l 2 a /(2b/ p 0 ) = t'° bs . At this time scale 

the initial shock front in the laser-generated acoustic strain pulse is well-defined and even slow 
shock front motion in the sonic coordinate frame could be still detected as it is described by Eq. 
(27) and Eq. (28), if the measurement of co dRIR are possible with the necessary precision. Note, 

that (co dRIR - co B ) I co B oc sM a (t) oc v sh (t) I c 0 could be very small in this quasi-linear regime. 

Consequently, the monitoring of the shock front motion could be possible due to the creation of 
weak shock front by pump laser action and not in the process of the acoustic pulse propagation. 
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The above presented theoretical analysis demonstrates that through the measurements of 
the frequency co dRIR of the time -resolved reflectivity signal there is an access to the measurements 
of the parameter of acoustic nonlinearity, because the nonlinearity is responsible for the difference 
between co dR/R and co B . Speaking differently co dRIR - co B 0 only if s ^ 0 . Although the 

information on s could be formally extracted from the experimental data in any regime of the 
nonlinear propagation of the coherent acoustic pulse, the most straightforward analysis of the 
experimental data can be achieved in the regime s Re » 1 when in the time interval 
2(2kl a y 2 t sh <t<2(sRd4) 2 t sh the transient reflectivity oscillation is caused by light scattering 
from a single well-defined moving shock front in the profile of the coherent acoustic pulse. In this 
regime co dRIR - co B = 2kv sh {t) . For 2&/ fl (0)«l in this regime the relations t»t sh and 

■Jl + (t -t sh )/(2t sh ) ~ yjt 1 2t sh hold, resulting in a particularly simple prediction for the time 

evolution of the measured frequency co dRIR -<x> B = \ — — co B <x -JsS It ccl/Jt. Consequently, 

V 2( V 

through the fitting of the experimentally measured [co dR/R (t)-a) B ]/ co B only the product 
so a x a = so a l a I c Q =sS can be measured. To extract the nonlinear acoustic parameter it is 
additionally necessary either to measure o a z a in some other experiments or to estimate o a z a 
theoretically. 

The dependence of the predicted co dRIR - co B in the considered regime only on the initial 
surface area of the acoustic pulse S = o a l a = o a r a c 0 , but not on o a and z a separately, indicates that 

the information on the initial pulse shape has disappeared at the considered long time delays. From 
this point of view fitting of the experimental data at shorter time scales, although more 
complicated, could be more fruitful. For example, in the theoretical prediction in Eq. (27) the 
variations of phase in time depend not only on the displacement of the shock front % sh {t) but also 
on the initial pulse length (duration) through the dependence on the non-dimensional 
parameter A(t) = 8^fc[/ a + g sh (t)/2] = 87ik[c 0 T a + <^ sh (t)/2] (see Fig. 3). Correspondingly, 

co dRIR -co B in Eq. (28) depends on the weak, i.e., acoustic, shock front velocity 



i sot 

1 + — — 



and K{t) = %M a (t) = &7ikc 0 T o 



so J 

1 + 



, and, thus, it is 



the function of two independent parameters so a and z a . Both parameters can be determined by 
fitting the experimental data. Then to extract the parameter of the acoustic nonlinearity it is 
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necessary to determine through additional experiments or estimate theoretically only the amplitude 
v a of the particle velocity in the photo-generated acoustic pulse. This situation is similar to one 

existing in the physics of strong, i.e., non-acoustic, shock waves, where in order to determine the 
equation of state at high pressures in addition to the measurements of the strong shock front 
velocity it is necessary to measure independently the flow velocity, i.e., the velocity of the particles, 
at least at one particular coordinate or in one particular moment of time. 31 " 36 

3.5. Brillouin scattering caused by multi-phase nonlinear acoustic pulses in the regime 
of arbitrary Reynolds numbers 

The most complicated for the approximate analytical description is the propagation of the 
acoustic pulses characterized by the condition e Re oc 1 , when the roles of the nonlinear and the 
dissipation effects in the Burgers equation (Eq. (10)) are comparable from the beginning. In this 
case there is in principle formal opportunity to describe the evolution of the pulse profile precisely, 
however thus obtainable solutions in most of the cases are too cumbersome to provide insight in 
the physics of the processes of interest 23 . To obtain the analytical solution of Eq. (10) it is first 
convenient to rewrite it by introducing non-dimensional slow time 9 = t I 't sh = so J I l a , non- 
dimensional profile coordinate a = %ll a , non-dimensional particle velocity V = v I v a and the 
compact notation l/(^Re) = T/4. The equation takes the form 



which confirms the crucial role played by the parameter sRc, as it had been discussed 
earlier. The nonlinear transformation of the function, the so-called Hopf-Cole transformation, 23 



dV_ 



+v 



dV r d 2 V 
da ~Ao~a i 




(33) 



V = - 



r dint/ 
2 da 



r\ l du 

2JU~^a , 



(34) 



transforms the nonlinear Eq. (32) in linear diffusion equation 
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dU r d 2 U 



d& 4 da 2 



= 0. 



(35) 



The diffusion equation has formal analytical solutions in the integral form, for arbitrary initial 
conditions U(0,a) 



,\2 



co f a-CT 
a J 



U(3, a) = -j— [e [ ^ ) U(0, a')da' = [e 1 a } U(0,a')da 



(36) 



Here we have introduced the compact notation a{&) = vr>9 . The solution for the spatial derivative 
can be presented, using Eq. (34), in the convenient form 



,\2 



cU(9.a) 1 J^J dU(0,a') dy _ -1 



da 4na 



da' 



-Jna 



j\2 



e K a > U(0,a')V(0,a')da'. (37) 



The initial distribution U(0,a) can be obtained from the initial conditions V(0,a) by using the 
inversion of Eq. (34), i.e., 



U{9, a) = exp 



r 



\v{9,a')da' 



= C exp 



v 

-^V($,a')da' 







= C exp 


^S(&,a) 





(38) 



Here C is the integration constant and the function S(a) has the sense of the surface area of the 
acoustic pulse starting from the co-ordinate a up to the pulse beginning at + go , i.e., to the right of 
the co-ordinate a . The more independent parameters are necessary for the description of the initial 
acoustic pulse profile V(0, a) , the less are chances to get the solution of Eq. (35) not in the form 
of the integrals or infinite sums, but in functional form on the basis of common and not "exotic" 
special functions. 

An important in picosecond laser ultrasonics class of the initial strain profiles corresponds to 
the generation of the acoustic pulses in the thin light absorbing films, when the laser-induced 
heating of the film material can be considered as spatially homogeneous and instantaneous, while 
the cooling can be considered as infinitely long. What is emitted from the film in the surrounding 
material is the sequence of the rectangular strain pulses with their amplitudes and signs controlled 
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by the reflection coefficients of acoustic waves at the interfaces between the film and the 
surrounding material. The duration of each pulse in the sequence is equal to the time of sound 
propagation across the film x a = H f I c f , where H f is the film thickness and c f is the sound 

velocity in the film. The typical particle velocity profile launched outside the film is presented 
qualitatively in Fig. 4. Although the positions of the successive fronts (i=l, 2, 3, ...) in the 
initial profile in Fig. 4 are equidistant, separated by the same distance l a =cr a ={cl c f )H f , we 

have not included this fact in the notations in Fig. 4 for the generality. When the normalized 
coordinate a = %ll a , the normalized particle velocity V = v I v a and the normalized time by 

9 = tlt sh =tl[l a l{eu a )\ are introduced, the normalized profile in Fig. 4 provides the initial 

condition for Eq. (33). Then the initial condition for Eq. (35) can be obtained, using Eq. (38), in 
the form 



U(0, a n+l < a < a n ) = U n (0, a) = C exp {(2 / T)[S n + V n (a n -a)]}, 



(39) 



where S n = S(0,a n ) = ^V j _ 1 (a j _ l - <jj) is the surface area of the pulse between a = a n and 



7=0 



a = co = a 0 and V(cr = 00) = ^ =0. With this boundary condition, using Eq. (36) and Eq. (37), 
we derive 



1 °° 

U(3,a) = ^U n (0,a)e^ 2 \erf 

2 „=o 



a„ -a + aV/r 



a 



erf 



a n+l -a + a%/r 
a 



dU(3,cr) 
da 



1 



ZV n U n (0,a)e^ 2 lerf 



n=0 



a -a + a 2 VJT 



a 



erf 



(40) 



Here U n (0,a) is defined by Eq. (39). Substitution of Eq. (40) into Eq. (34) provides analytical 

description of the acoustic propagation for arbitrary magnitude of the parameter ^Re = 4/r. 
However the Fourier transform of the predicted signal, when it is substituted in Eq. (11) for the 
evaluation of the oscillations in transient optical reflectivity, has to be done numerically. As a 
simplest example we present the solution in Eq. (40) for the initial single rectangular pulse 
localized between a 2 =-1/2 and cj x =1/2 and characterized by a single non-zero amplitude 
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U{9,cr) = ^\\-erf 



l/2-o- 



+ 2<\I2-<j)IT (aVJT)' 



-erf 



f -l/2-a + a 2 VJT 



y a J 



erf 



l/2-a + a 2 VJT 



J 



+ e 



2/r 



dU (>9, cr) _ e 2(i/2-cr)/r g (av„ irf 
da ~ T 



■1/2-er 



+ 1 



(41) 



^1/2-a + aV/r^ 



' -l/2-a + a 2 K/r^ 



Using either the solution in Eq. (41) or general formulas in Eq. (36) and Eq. (37) in the case 
of the initially symmetrical pulse, it is straightforward do verify that the asymmetry in the pulse 
profile, which is introduced by the nonlinearity when the pulse starts to propagate, inevitably 
disappears at longer time scales. For example, the general solution in Eq. (37) for the derivative 
takes, in the case of the initial rectangular pulse at the time scales a{&) > 1 , the following form 



dU(3,cr) 



-Jna 



1/2 o"- 0 " 
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-1/2 



U(Q,a')dcr'z 



■ina 



1/2 



-1/2 



+ 2 
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.rl 2 



} d&. 



The integration leads to the following ratio P between the asymmetric, proportional to a/ a, and 
symmetric, independent of the co-ordinate, contributions to the derivative: 
P * -[tanh 1 (1 / r) - r\a I a{3)]. Thus in the limit of strong absorption, when ^Re«l, while 
r»l, we estimate P « -(l/3)(l/r)[cr/a(<9)] . This indicates the low level of asymmetry as 
expected. However even in the regime of the initially strong nonlinearity when s Re »1, while 
r«l, we estimate P « -[cr/a(<9)] . This indicates the diminishing of pulse asymmetry with 
increasing time when a{&) > 1 and that the asymptotic of the evaluated derivative will be the 
Gaussian profile. At the same time the asymptotic form of the profile U (9, a) , which follows from 
Eq. (41) when a(3) »1, is U(3,&) = C[l + exp(2/r)], that is co-ordinate independent. Thus in 
accordance with Eq. (34) the asymptotic form of the profile of the acoustic pulse V(9 — > oo, cr) is a 
Gaussian, confirming the statement formulated in Section 5.1 on the disappearance of the acoustic 
pulse asymmetry at sufficiently long time scales . 

The importance of the solution in Eqs. (34) and (40) is in the opportunity, by applying in Eq. 
(11) the numerical Fourier transform, to derive the predictions for the evolution of the amplitude 
and the phase of the Brillouin oscillation in the most complicated regime characterized by 
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s Re cc 1 , where the asymptotic methods could be not enough effective and precise. The 
experimental data for co dRIR -co B should be fitted by dq> n I dt , which follows from Eqs. (11), (34) 

and (40), by using the fitting parameter sRe. When the initial pulse profile is well described 
theoretically while the effective viscosity b is known from the measurements of the linear sound 
absorption and velocity, the fits provide the value of so a . Finally, to estimate the nonlinear 

parameter s , the acoustic pulse amplitude v a should be measured independently or estimated 
theoretically with the knowledge of the linear physical properties of the materials. 

3.6. Brillouin scattering by the bipolar coherent acoustic pulses 

The solution presented in Eqs. (11), (34) and (40) is valid for the sequences of the 
rectangular acoustic pulses of arbitrary polarity. The particular solution in Eq. (41) describes the 
evolution of a single uni-polar acoustic pulse, which can be in practice emitted from a thin light 
absorbing film in the cases when the film is either perfectly matched in terms of acoustic 
impedance with the surroundings at its both faces, or when the film is perfectly matched with the 
substrate on one face while the other face of the film is loaded by a material with much higher 
(infinite) impedance. In the limiting case ^Re »1 the evolution predicted by Eqs. (34), (41) and 
the Brillouin signal predicted by Eq. (11) are qualitatively similar to those predicted for the uni- 
polar acoustic pulses by the asymptotic methods of the analysis in Sections 3.1 -3.4. 

Another common experimental situation is the case where one face of the film is in 
contact with the substrate while the other face is mechanically free, i.e., it is loaded by air with 
negligibly small acoustic impedance. In the limiting case of the perfectly matched film and the 
substrate, the coherent acoustic pulse launched in the substrate is bi-polar composed of two 
rectangular pulses of the opposite polarity. The analysis of the evolution of this pulse in the regime 
^Re»l by the asymptotic methods described in Sections 3.1-3.4 provides insight on what 
additional phenomena could be expected, at least qualitatively, in the case of an arbitrary sequence 
of the rectangular pulses of varying polarity described in Section 3.5. We will not present here the 
complete approximate analysis which can be conducted along the lines presented in Sections 3.1 - 
3.4, but we will just mention the most important difference between the uni-polar pulses (see Fig. 
2) and the multi-polar acoustic pulses (see Fig. 4). 

In the multi -polar pulses multiple shock fronts could be developed in the regime s Re » 1 or 
just exist initially due to the specific opto-acoustic generation process. Then it could be expected 
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that the time-resolved reflectivity measurements will be sensitive to all shock fronts and, generally, 
the time-resolved Brillouin oscillation will be more complicated due to the interference between 
the contributions to the total signal from all shock fronts. This is true except for perfectly 
symmetric bi-polar acoustic pulse, which presents opportunity for drastic simplification of the 
results due to the assumed symmetry. For example, at times scales, when for a uni-polar pulse in 
Sections 3.1-3.4 the Brillouin signal in Eq. (25) is due to probe scattering by a shock front with the 
amplitude su a (t)l2 located in the sonic coordinate frame at £ = £ sh (t), while the light scattering 
from the trailing edge of the pulse is negligible, the signal from the bi-polar pulse will contain 
addition contribution accounting for the second shock front with the amplitude -so a {t)l2 

located at £ = -% sh it) . The spectrum of the strain can be obtained using Eq. (26) with 
do I d{ = -v a [S{Z + £ sh (0) + S(4 + U (f ))] 

r/ zz (t,q = 2k) = i^cos[2*£ A (0]. (42) 

KCq 

The solution in Eq. (42) reveals an important qualitative difference which could be expected 
in some cases of multi-polar pulses in comparison with the case of a uni-polar acoustic pulse. It 
can be seen that in the case of a perfectly symmetric bi-polar pulse the nonlinear processes and the 
shock front motion do not manifest themselves by the variations in frequency of the reflectivity 
oscillations because the unrapped phase of the signal (42) does not depend on time, <p = const , 

and, consequently, co dRIR = co B . The nonlinear processes manifest themselves as an additional 

modulation, proportional to |cos[2fc£ sA (?)] , of the Brillouin amplitude, which is superimposed on 

the asymptotic dependence \(t) <x v a {t) . Thus the access to the nonlinear parameter in this case 

is through the measurements of the variations with time of the amplitude of the oscillation in 
transient reflectivity signal. It is worth mentioning that in the regime, currently under consideration, 
2k^ sh (t) coincides with the parameter A(0 introduced in Eqs. (27) and (28) in Section 3.2. In 

accordance with Eqs. (27) and (28) the additional modulations of the reflectivity signal due to the 
light scattering by the trailing front of the positive phase and the leading front of the negative 
phase of a bi-polar pulse can be important when A(t) < 1 . For A(t) > 1 the oscillatory character of 
the amplitude modulation in Eq. (42) could be important. So for the convenience it is worth to 
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combine these predictions in a single formula. For the bipolar pulse described 
v = v a (f )[(£ / U (t)]H(Z + £ A (f))[l - - l h (t))] we derived from Eq. (1 1) that 

1^1 = 2k) = ^{cosfeCO]-*^}. (43) 
kc 0 { 2^(0 J 



4. Comparison of the theoretical predictions with the experimental observations 

In this Section we apply the developed theory to the interpretation of results obtained 
from conventional picosecond laser ultrasonics experiments. In these experiments the formation 
of weak shock fronts was expected during the propagation of the picosecond acoustic pulses. 

4.1. Spatially resolved Brillouin scattering technique 

The time-dependent amplitude oscillation of a time-resolved Brillouin scattering signal as 
predicted by Eq. (43) for bi-polar nonlinear acoustic pulses can be measured experimentally 
though observing the oscillation in space. This is possible because the modulation of the Brillouin 
scattering amplitude in time is equivalent to its modulation in space. Therefore, it is sufficient to 
substitute r => z/c 0 in the derived theoretical formulas. Recently, these spatial modulations have 

been experimentally observed in picosecond laser ultrasonics by detection of ordinary, i.e., 
stationary, Brillouin scattering of continuous light at multiple spatial positions along the 

24 

propagation path of coherent picosecond acoustic pulses propagating in a sapphire sample. In 
these experiments, a coherent strain pulses was generated in a 100 nm chromium film with one 
mechanically free surface. For such a 100 nm chromium film the characteristic initial duration of 
the acoustic pulse is r a = H f I c f « 17 ps. The nonlinearity of sapphire is characterized by the 

following nonlinearity constant 37 : 2p 0 c%£ = 18.3xl0 n N/m 2 (f»1.82); while the effective 

viscosity at room temperature is assumed to be 24 & = 4.54xl0 4 sN/m 2 . For room temperature 

24 —3 

experiments and using a strain amplitude of M a = v a I c Q « 10 we estimate ^Re«68»l. 

From this estimate we conclude that the results of this experiment can be analyzed by the 
approximate methods described in Section 3. The characteristic nonlinear time is 
t sh = r a l{eM a ) « 9.3 ns and the time it takes for the initially rectangular compression phase of the 
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photo-generated pulse to transform into triangular form is 2t sh =~ 18.6ns. It could be expected that 
the two shock fronts will be well defined at the times significantly shorter than 
= 2(sRd 4) 2 t sh « 5S0t sh « 5.4 //s. Thus, we expect the amplitude modulation as described by 
Eq. (43) for the time window 20 ns < t < 500 ns . Now, when we consider the sound velocity in 
sapphire as c 0 = 11.23xl0 3 m/s, this time window will correspond to pulse propagation distances 
from approximately 0.2 mm to 5 mm, which exactly corresponds to the interval studied 

24 

experimentally by Muskens and Dijkhuis . For the used probe wavelength of 514 nm, the 
amplitude of the transient reflectivity signal starts to drop faster than predicted by Eq. (43) at the 
end of this observation interval because of the broadening of the shock fronts as described by Eq. 
(31). However, both, our analytical estimates and the numerical modeling of the solutions of the 
Burger's equation (considering the conditions from Ref. 24), confirm the existence of well-defined 
fronts for bipolar acoustic pulses up to a propagation distance of 5 mm. Thus, we attribute this 
experimentally observed oscillatory modulation of the scattered-light amplitude to the fact that 
under these conditions the light is predominantly scattered by two shock fronts moving in the sonic 
co-ordinate frame with equal, but oppositely directed velocities. The oscillatory modulation is a 
result of the continuous increase in the phase difference between these two signal contributions. 

4.2. Time-resolved Brillouin scattering technique 

More recently 25 , results from time-resolved Brillouin scattering of nonlinear picosecond 
acoustic pulses in liquid were reported. Note that this is the detection technique explicitly assumed 
in our preceding theoretical analysis in Sections 2 - 6.. In those measurements, the photoexcitation 
of a H=33 nm thick aluminum film deposited onto a sapphire substrate and loaded by a liquid, 
DC704, ' launched a train of acoustic pulses of alternating polarity and progressively 
diminishing amplitude into the liquid. The structure of the pulse is illustrated in Fig. 4. The 
duration of each of the alternating phases is equal to twice the time of longitudinal acoustic wave 
propagation across Al film. Thus, the characteristic duration of the phases at half-width-half- 
maximum can be estimated as r a = H I c Al ~ 5.2 ps, where c Al « 6.4 -10 3 m/s is the sound velocity 
in Al. The relative amplitudes of the consecutive parts of the emitted acoustic wave are controlled 
by the reflection coefficients of sound at the Al-sapphire and Al-liquid interfaces , R { = R A u sapphire 

and R 2 = R AUDC1(H , respectively. Using R l =0.45 , R 2 =-0.75 , we estimate the shape of the 
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acoustic pulse shown in Fig. as v a2 1 v al = R x = 0.45 , o a3 1 v al = R X R 2 = -0.34 , 
v aA I v aX = R?R 2 s -0.15 , v a5 1 u al = R?R? =0.12. 

The characteristic time t d at which the thermal diffusion in Al is supersonic 2 ' 28 can be 
estimated as t d = D M I C 2 M « 1.1 ps, where D A[ = 4.5 ■ 10 5 m 2 /s is the thermal diffusivity of Al. At 
times t > r D the thermal diffusivity is subsonic and has weak influence on the profiles of thermo- 

1 2 28 25 

elastically generated acoustic pulses. ' ' Under the conditions of the experiment r D is an order 

of magnitude shorter than the full duration, 2r a =2Hlc M «10.4ps, of each of the alternating 

acoustic phases. That is why it is reasonable in the first approximation to neglect the contribution 
to the leading compression phase of the emitted pulse train in Fig. 4 of the acoustic waves emitted 
in the process of the Al film cooling, as it had been apriori assumed when plotting Fig. 4. The 
signal profile in Fig. 4 includes only the acoustic waves generated due to laser heating of Al film. 

In the experiments the probe laser radiation had been incident on Al film from the liquid 
side and the detected oscillating in time contribution to probe reflectivity dR(t) I R is mostly due 
to the interference between the probe radiation reflected by the Al film and reflected by the 
hypersound launched in liquid. The Brillouin frequency f B = co B I{2tt) = 40 GHz had been 
measured at low pump laser fluence F L < 0.2 mJ/cm , where it had been observed that the 
frequency co dR/R of the oscillating component of transient reflectivity does not depend on pump 
laser fluence. The major result of the experiment 25 is the reported dependence of the frequency 
to dRiR °f me oscillating component of the transient reflectivity dR(t)IR both on pump laser 
fluence and on time in the experimentally tested ranges 0.29 mJ/cm2< F L < 10.3 mJ/cm 2 , 100 ps 
< t < 2000ps, corresponding to more than 35 times increase in laser fluence and an order of 
magnitude wide window in the time domain. 

In accordance with the theory developed above in Section 3 the analysis of the experimental 
data should be started by estimating the characteristic parameter s Re , the product of the nonlinear 
acoustic parameter and of the acoustic Reynolds number, which determines the regime of the 
emitted acoustic wave evolution in time and space. In the experiment the relative deviation 
(to dR/R -co B )l co B of the detected oscillation frequency from the linear Brillouin frequency in the 
beginning of the evaluated time interval, i.e., at f=100 ps, is reported to be 
0.002 < {co dRIR -co B )l co B < 0.072. These values can be used to estimate the product sM a of the 
acoustic nonlinear parameter and of the acoustic Mach number, by assuming that the deviation of 
to dR i R from co B are caused by the supersonic propagation of the weak compression shock related to 
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the leading compression phase of the acoustic signal presented in Fig. 4. In accordance with Eq. 
(25), we use for the estimates the relation (co dR/R -co B )lco B »sM a (t)/2. Consequently, at t=!00 

ps we can estimate 0.004 < sM a < 0.144. Using the definition of the acoustic Reynolds number, 23 

Re = 2p Q c 0 u a r a I b , and the knowledge of sound absorption coefficient a = a(a>) in DC704 at 

25 2 3 1 

considered frequencies and temperatures a(co)=bco /(2p Q c Q ) =T I c Q , with T = 2 ns" at the 
experimental temperature of 200°K, we present sRe in the form sRe = sM a co 2 B T a IT . Thus for 
0.29 mJ/cm 2 <F L < 10.3 mJ/cm 2 at t=100 ps it can be estimated that 0.63<^Re<23 . It is 
important that these estimates demonstrate that for F L > 4.4 mJ/cm the non-dimensional 
parameter <?Re> 10»1 and, theoretically, the evolution of the leading compression phase of the 
acoustic signal in Eq. (4) can be described not by the Burgers equation, Eq. (10), but by the 
simple-wave equation, Eq. (14), which describes the dynamics of delta-localized shock fronts. 

In accordance with the theory presented in Section 3 the largest in amplitude weak 
compression shock front is expected in the evolution of the leading compression phase, for which 
the above estimates have been conducted. The largest in amplitude dilatation shock front could be 
formed at the trailing edge of the second (dilatation) phase of the acoustic signal in Fig. 4. 
However, in accordance with our estimates the amplitude of the dilatation shock would be at least 
° a i /(°a5 ~ v aA) ~ 4 times smaller than of the compression shock, the condition s Re»l does not 
hold even at the highest pump laser fluences, and, consequently, the dilatation shock will be 
importantly additionally suppressed relative to the compression one by the acoustic absorption. 

The conducted estimates indicate that under the conditions of the considered experiment it is 
worth trying to model transient Brillouin reflectivity signal by assuming that it is due only to the 
first leading compression phase of the acoustic pulse train. Moreover, assuming that the 
compression pulse in the interval £ 3 < £ < 0 (see Fig. 4), which is composed of two rectangles, is 

transformed due to nonlinear acoustic process in a triangular shape pulse at time scales r sh < 100 

ps, the experimental data can be fitted in the complete observation window 100 ps < t < 2000ps by 
the formulas derived in Sections 3.1 - 3.4 for the triangular pulses. Neglecting in Eq. (28) the 
modulation of (co dR/R -co B )l co B due to the probe light reflection from the trailing edge of the 
triangular pulse, the fitting formula can be presented in the form 
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Here M a =v a l c 0 is the initial, i.e., at t=0, Mach number of the equivalent triangular pulse. 
The initial amplitude v a of the triangular pulse is related to the initial amplitude v al of the leading 
phase in the signal presented in Fig. 4 by the law of the pulse area conservation, i.e., 
v a i T a + v a2? a = (1 + ^ i )^ a i T a = v J a ■ Consequently, v a = (1 + R x )u al and 

sM a =(\ + R l )sM aX K\A5sM aX . The fits with Eq. (44) were conducted [25], by adjusting the 
parameters sM a and r a , all the experimental dependences measured for 0.29 mJ/cm 2 < F L < 10.3 

mJ/cm . Even visually the quality of the fits continuously increases with the increasing pump laser 
fluence due to the increasing characteristic acoustic Reynolds number, as it could have been 
previewed from the theory and the estimates. At fluences F L > A A mJ/cm , where validity of Eq. 
(44) is expected, the dependence of the acoustic Mach number and the initial duration of the 
compression pulse can be well approximated by 

sM a * 0.035(F L /[mJ/cm 2 ]) , zf = 5.5 ps . (45) 

The close proximity of the rf' to the theoretically expected value z a =H I c Al ~ 5.2 ps is a strong 

support for the validity of the theoretical model in Eq. (44). Moreover, using Eq. (45) it is 
straightforward to estimate that the characteristic time of the triangular pulse formation, 
r Sh = (c 0 T a )/(su al ) = (1 + R l )r a l(sM a ) , is shorter than 100 ps for F L > 2.4 mJ/cm 2 . It is also 

possible to estimate, by using Eq. (30), that the finite width of the shock front, A^ sh (t) , which is 
controlled by the competition of the nonlinearity and dissipation, is shorter than the initial duration 
r a of the acoustic pulse even at the end of the experimental time window, i.e., at t=2000 ps, for 

F L > 4.3 mJ/cm . Finally, it should not be forgotten, that in the model of Eq. (45), in comparison 
with more precise description in Eq. (28), the modulation provided by the function S(A(0) and 
presented in Fig. 3 (b) has been omitted. In accordance with the plot in Fig. 3. (b) the effect of this 
additional modulation is smaller than 10% for A(t) > n . This inequality is equivalent to 

T a(0 = ? a -\A + ( £M a 1 2 )( ? I T a) ^ 7i l{2o> B ) ~ 6.2 ps. Using Eq. (45) we estimated that the additional 
modulation provided by the probe light reflection from the trailing edge of the compression pulse 
is negligible at time scales t > 100 ps for pump fluences F L > 1.9 mJ/cm 2 . 
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On the basis of the above presented estimates and good quality of the fits it can be concluded that 
the deviation of the frequency co dRIR of the oscillating contribution to transient reflectivity from the 

linear Brillouin frequency co B for F L > 4.4 mJ/cm in the experiments [25] is due to the 
supersonic propagation of the weak compression shock front. Moreover the dynamics of the weak 
acoustic shock front evolution at picoseconds - nanosecond time scale can be described by the 
classical formulas of the nonlinear acoustics. The experimentally determined dependence in Eq. 
(45) of the product sM a on the pump laser fluence can be used for the evaluation of the nonlinear 

acoustic parameter of liquids in GHz frequency range if the dependence of the acoustic Mach 
number on laser fluence is modeled theoretically from the thermo-elastic sound generation 
problem or measured experimentally (for example by monitoring optically the displacement of the 
Al/sapphire interface at picosecond time scales). 

5. Conclusions 

We presented an analytical theory describing how the nonlinear reshaping of a propagating, 
finite amplitude picosecond acoustic pulse results in the time-dependence of the observed Brillouin 
frequency and amplitude in a time-resolved Brillouin scattering experiment. In particular, the 
conditions, under which the detected frequency provides information on time-evolving speed of 
the weak-shock front or fronts that can be formed through the nonlinear acoustic transformation of 
the strain pulse, are established. Depending on the parameters of the initial strain pulse and the 
monitored time interval of its nonlinear transformation, our theory predicts the detected frequency 
to either be monotonically decreasing or oscillating in time. We also predicted the influence on the 
Brillouin scattering signal of the linear acoustic absorption which could broaden with time/distance 
the widths of the propagating shock fronts. We highlighted different behaviours of the Brillouin 
signals in case of uni-polar and bi-polar nonlinear acoustic strain pulses. We presented examples 
how the developed theory can be successfully applied for interpretation of available experimental 
data, providing access to the evaluation of the nonlinear acoustic properties of the media in the 
GHz frequency range. 
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Figures and captions 




Figure 1 

Physical principle of the picosecond acoustic interfermetry. Weak optical field reflected by the 
plane acoustic strain pulse propagating normally to the sample surface is heterodyned against 
stronger field reflected from the sample surface. When the acoustic pulse propagates at constant 
velocity, linearly changing in time phase shift between the two interfering reflected optical fields 
leads to sinusoidal variation with time of transient reflectivity changes. 
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Figure 2. Transformation of the particle velocity profile of the triangular acoustic pulse of the 
finite amplitude in the frame of the quadratic approximation of the nonlinear acoustics, (a) All the 
points of the initial profile with the peak located in point A are displaced with increasing 
propagation time horizontally. These displacements are illustrated for some particular points by the 
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arrows. This leads first to the shortening of the leading front of the acoustic pulse (profile with the 
peak in point A' ) and later to formation of the shock front (profile with the peak in point A"), (b) 
After the formation of the shock front (profile with the peak in point A" ) the amplitude of the 
acoustic pulse amplitude starts to diminish with increasing time (profile with the peak in point A'" ). 
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Figure 3 

(a) The presentation of the function H(r) = *Jl - 2 sin c(A) + sin c 2 (A / 2) , describing transitional 

modulation of the Brillouin signal amplitude, which is superimposed on its asymptotic 
proportionality to the velocity of the shock front. 
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(b) The presentation of the function 0(A) = smc_ ^ describing transitional 

l-2sinc(A) + sinc (A/2) 

modulation of the frequency co dRIR of the transient reflectivity oscillation, which is 

superimposed on its asymptotic variation proportionally to the velocity of the shock front. 

Here A(t) = 4kl a (t) is the normalized length of the acoustic pulse, which starts to increase after the 

shock front formation. The modulations are due to the contribution of probe light scattered by the 
trailing long front of the acoustic pulse. 
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Figure 4 

Presentation of the typical train of rectangular acoustic pulses emitted into the sample from thin 
metallic film sandwiched between the substrate and transparent loading material and heated by 
ultrasort laser pulse. It is assumed, for definiteness, that the acoustic impedance of the material is 
lower than that of the film, while the impedance of the film is lower than that of the loading. 



